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ABSTRACT 



Context. The observations carried out by the space missions CoRoT and Kepler provide a large set of asteroseismic data. 
Their analysis requires an efficient procedure first to determine if a star reliably shows solar-like oscillations, second to 
measure the so-called large separation, third to estimate the asteroseismic information that can be retrieved from the 
Fourier spectrum. 

Aims. In this paper we develop a procedure based on the autocorrelation of the seismic Fourier spectrum that is capable 
of providing measurements of the large and small frequency separations. The performance of the autocorrelation method 
needs to be assessed and quantified. We therefore searched for criteria able to predict the output that one can expect 
from the analysis by autocorrelation of a seismic time series. 

Methods. First, the autocorrelation is properly scaled to take into account the contribution of white noise. Then we 
use the null hypothesis Ho test to assess the reliability of the autocorrelation analysis. Calculations based on solar and 
CoRoT time series are performed to quantify the performance as a function of the amplitude of the autocorrelation 
signal. 

Results. We obtain an empirical relation for the performance of the autocorrelation method. We show that the precision 
of the method increases with the observation length, and with the mean seismic amplitude-to-background ratio of 
the pressure modes to the power 1.5±0.05. We propose an automated determination of the large separation, whose 
reliability is quantified by the Ho test. We apply this method to analyze red giants observed by CoRoT. We estimate 
the expected performance for photometric time series of the Kepler mission. We demonstrate that the method makes 
it possible to distinguish £ = from £ — 1 modes. 

Conclusions. The envelope autocorrelation function (EACF) has proven to be very powerful for the determination of 
the large separation in noisy asteroseismic data, since it enables us to quantify the precision of the performance of 
different measurements: mean large separation, variation of the large separation with frequency, small separation and 
degree identification. 
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1. Introduction 

Astcroscismology is known to be an efficient tool to ana- 
lyze the stellar interior and to derive the physical laws that 
govern stellar structure and evolution. It benefits nowa- 
days from high-performance photometric data provided by 
the space missions CoRoT ( Baglin et al. 200"6| ) and Kepler 
( jChristensen-Dalsgaard et al. 2007 ). The amount of data is 
much higher than from the earlier ground-based observa- 
tions, even with the recent multi-site ground-based observa- 
tions ([Arentoft et al. 2008|) . since space-borne instruments 
are able to simultaneously record long time series on nu- 
merous targets. The data analysis then must be efficient 
enough to rapidly extract seismic information from hun- 
dreds to thousands of stars. 



* The CoRoT space mission, launched on 2006 December 27, 
was developed and is operated by the CNES, with participation 
of the Science Programs of ESA, ESA's RSSD, Austria, Belgium, 



This task is principally carried out on the frequency pat- 
tern of the eigcnmodes propagating inside the stars. For 
targets showing solar-like oscillations, this pattern follows 
the asymptotic relation of Tassoul (1980) providing eigen- 
frequencies nearly equally spaced by Av/2. The eigenfre- 
quency of radial order n and degree £ expresses v n> i ~ 
[n + £/2 + s]Av - £(£ + 1)D , Av being called the 'large 
separation, Dq giving a measure of the small separation, 
and e a constant term. The determination of the large sep- 
aration Av is the first step of any seismic analysis. If the 
signal-to-noise ratio is high enough, Av can be detected by 
eye in the power spectrum. In many cases, this is not pos- 
sible, and the determination of Av requires sophisticated 
tools, as was the case for the first correct determination 
of the large separation of Procyon (jMosser et al. 1998)) and 
of the first CoRoT target observed with Doppler measure- 
ments ([Mosser et al. 2005|) . For observations dealing with 
a single target, the tools used for the determination of Av 
are usually unautomated and involve parameters specific 
to the target. Most often, they require the visual inspec- 
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ting in other cases. Discussion of various cases is presented 
in Section We propose an automated determination of 
the large separation; using the Ho test, we can quantify the 
reliability of this method. Section [6] is devoted to conclu- 
sions. 



tion of an image or a graph obtained by transforming the 
Fourier spectrum using the asymptotic relation cited above 
(echelle diagram or comb response) . This step can be auto- 
mated, but with great care, since the higher order terms of 
Tassoul (1980) complicate the stacking, as does, for exam- 
ple, the variation of the large separation reported in many 
asteroseismic targets (e.g. IMosser et al. 2008[) . 
With the advent of space photometric missions, the use of 
pipelines for the automatic detection of Ai> is becoming 
mandatory (IMathur et al. 20"09|) since many targets have 
a low signal-to-noise ratio. A test to determine if the large 
separation is reliably detected is highly desirable, and a way 
to estimate the asteroseismic content of a high-precision 
photometric time series will be very helpful. 
This paper proposes an original way to address 
these issues. It is based on a first report by 

(hereafter RV06), who C(t) = x(t) x(t + r)dt = X(v)X*(v)e 



2. Autocorrelation 

2.1. Calculation 

RV06 proposes to perform the autocorrelation of the seismic 
time series as the Fourier spectrum of the filtered Fourier 
transform of the time series. This directly gives the am- 
plitude of the envelope of the autocorrelation function, as 
shown in the Appendix. Instead of the canonical form, 



Roxburgh fc Vorontsov (2006) 
analyse solar-like oscillations via the square of the auto 
correlation of the time series, calculated as the Fourier 
spectrum of the filtered Fourier spectrum. RV06 state that 
the method is useful when faced with low signal-to-noisc 
ratio data, and might be useful in obtaining information 
about a s tar even when ind ividual frequencies cannot be 
extracted. Roxburgh (2009) (hereafter R09) shows that it 
is possible, with basic and rapid computations, to attain 
more complex objectives, such as measurement of the 
variations of the large separation with frequency. 
Since it provides a rapid measurement of the large sepa- 
ration, the autocorrelation method fits perfectly with the 
main asteroseismic objective of the Kepler mission, the 
large separation being used as an independent measure- 
ment in_e2ctractin£_theradius of stars hosting exoplanets, 
as in|Stello et al. (2009) The autocorrelation achieves this 



l2lTUT 



dp 



(1) 



with X{y) the Fourier transform of x(t), the autocorrelation 
with a filter T of width 8va centered on v c can be written: 

C= X{v)X*{v) T(v) c [2 ™ T di>. (2) 

We deal with the dimensionless square module of the auto- 
correlation: 



A* = \C{t) 2 \ / \C{0f 



(3) 



The choice of square module has no impact on the results 
presented below, but proved to be more convenient in many 
cases, such as the observed linear increase of A* with the 
observing time (see Eq. (JS])). 



goal without fitting a complex mode pattern to the stellar 2.2. Noise scaling 



power spectrum. Therefore, it provides a simple tool to es- 
timate the asteroseismic information of a Fourier spectrum 
or to use with Kepler, which will produce numerous time 
series of stellar targets. 

Wc propose to quantify the relevance of the autocorrela- 
tion method with the null hypothesis, and to determine 
simple criteria to assess its efficiency and predictive power 
when analyzing an oscillation spectrum with a low signal- 
to-noise ratio. The method is also useful for extrapolating 
the performance obtained with a short time series to that 
obtained with a 4-year long time series, as will be provided 
by the Kepler mission. The analysis relies on photometric 
time series as observed by CoRoT ( Baglin et al. 2006 1 , plus 
simulations based on these CoRoT spectra with the addi- 
tion of noise. It also includes simulations derived from a 
solar oscillation spectrum observed in photometry by the 
VIRGO/SPM instrument of the SOHO mission. 
Section [2] introduces the envelope autocorrelation function 
(EACF) and the way we scale it to properly account for the 
noise contribution. We show in Section [3] how the value of 
the main autocorrelation peak varies with different global 
parameters of the stellar oscillation spectrum. A crucial pa- 
rameter is the mean seismic height-to-background ratio 1Z, 
representing the smoothed height of the seismic power spec- 
tral density compared to the background. We introduce in 
Section 2] the Ho test, that allows us to examine and to 
quantify the performance of the method. The value of the 
EACF gives a reliable criterion to estimate the seismic out- 
put, from the determination of the mean large separation 



In order to compare different cases, it is preferable to ex- 
press the amplitude of the autocorrelation signal in noise 
units. The mean noise level in the autocorrelation can be 
derived from the fact that the noise statistic is a \ 2 with 2 
degrees of freedom. It is expressed in the general 



2_ <^> 



(4) 



with Nt the number of points in the time series. The noise 
level cr is inversely proportional to the number Na of fre- 
quency bins selected in the filter, when Na is measured 
in a Fourier spectrum at the exact frequency resolution 
dp = 1/T, T being the length of the observation, with- 
out oversampling. For a cosine filter (or Hanning function) 
of full-width at half-maximum Sva centered on v c , one gets: 



oh 



2iV H 



(5) 



With such a cosine filter and the resulting noise level, we 
define the EACF: 



A = A* /era. 



(6) 



We note that Aav = A{ta v ), the amplitude of the first 
peak in the autocorrelation function, at a time shift ta v = 

is 



2/Aza The first peak of the autocorrelation (Fig. IA.1[) 
the signature of the autocorrelation of a seismic wavepacket 
after crossing the stellar diameter twice. As shon in RV06, 
measuring the time delay tav of this peak allows us to 
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Table 1. TZ and parameters of the p mode envelope for CoRoT targets 
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HD49385 a 


GOIV 


7.41 


136.9 


1.00 
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0.95 


9.7 
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HD49933 61 


F5V 


5.78 


136.9 
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1.10 
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0.75 
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HD49933 62 






60.7 














237 


HD175726 C 


GOV 


6.72 


27.2 


2.05 


0.82 


97 


1.20 


8.5 


0.11 


6.9 


HD181420 d 


F2V 


6.57 


156.6 


1.60 


0.76 


76 


1.15 


10.3 


0.42 


242 


HD181906 e 


F8V 


7.6 


156.6 


1.92 


0.88 


85 


1.00 


10.3 


0.15 


47 


HD181907 / 


G8III 


5.8 


156.6 


0.0286 


0.0176 


3.5 


1.10 


5.0 


1.99 


40 


Sun/ VIRGO 9 


G2V 




182.1 


3.25 


1.04 


135 


1.30 


8.0 


2.50 


7.5 10 3 



i/ ma x is the location of the frequency of maximum power; Su cnv is the full-width at half-maximum of the mode envelope; (Av) 
is the mean value of the large separation; a represents the optimized filter width, in unit Svenv; 7 represents Svonv in unit Av; 
TZ measures the mean seismic amplitude in the time series compared to the noise, by the ratio in the Fourier spectrum, at the 

maximum-oscillat ion frequency, of the s mooth ed mode height to the b ackgr ound power density; yl ma x m easur es the EACF. 

IMosser et al. 20091 



References: tDehcuvcls ct al. 2009] ^Benomar et al. (2009)| ^Appourchaux et al. (2008) 
iBarban et al. 20091 IGarcia et al. 20091 ^Carrier et al. 20091 IFrohlich et all 19971 




Fig. 1. Contributions to the smoothed power density distri- 
bution for HD49933. The oscillation spectrum was slightly 
and severely smoothed (solid thin and thick lines). The 
dashed line represents the contributions of granulation and 
photon noise. The dash-dot lines account for the Gaussian 
modeling of the seismic envelope and the total contribution. 



3. Analysis 

We tested the variation of Aal> with various parame- 
ters, in order to determine the relevant ingredients con- 
tributing to this signal. We based the analysis on solar 
data obtained with the VIRGO/SPM instrument onboard 
SOHO (jFrohlich et al. 1997p . and on the CoRoT data pro- 
vided on the solar-like targets HD49933 (Appourchaux 
et al 2008), HD493 85 (jDeheuvels et al. 2009[), HD17 5726 
(jMosser et al. 2009|). HD18142 (jBarban et al. 2009|) and 
HD181906 (|Garcia et al. 2009|) . We al so include the red gi - 
ant HD181907 observed by CoRoT (jCarrier et al. 2009|) . 
All these targets are presented in Table [TJ We also consid- 
ered a set of red giants observe d in the exoplanetar y field 
of CoRoT, already analyzed by Hekker et al. (2009) 



3.1. Seismic amplitude-to-background ratio 

The strength of the autocorrelation of the time series de- 
pends on the ratio of the mean seismic amplitude compared 



noise ratio in the time series from the ratio estimated in 
the oscillation spectrum. This ratio in the Fourier spec- 
trum does not depend on the frequency resolution when 
the modes are resolved, i.e. the observation time is longer 
than the mode lifetime. In order to remove the influence 
of unknown parameters, such as the star inclination or the 
mode lifetime, we have to consider the ratio TZ of the mode 
height to the background power density, at the maximum- 
oscillation frequency, in a smoothed power density spec- 
trum (Fig. [!}. 

In order to estimate the background power, we have mod- 
eled the Fourier sp ectra with three components as in 
|Michel et al. (2008) a low-frequency Lorentzian-like pro- 
file, a Gaussian mode envelope and a high-frequency noise. 
Figure [1] shows this modeling for HD49933. The smoothed 
power density depends on the filter width. In order to avoid 
Gibbs phenomenon-like structures, a Gaussian filter has to 
be preferred to a boxcar average. The width has to be pro- 
portional to the large separation: a value of 3 Ai/ provides 
the optimum smoothing and limits the influence of the vary- 
ing background level. Since, at this stage, the large separa- 
tion is a priori unknown, the value of the filter width can 
be estimated with the help of the relation found between 
the large separation and the maximum-power frequency de- 
rived from the solar-like CoRoT targets: 



(0.24 ±0.05) v. 



0.78±0.045 

max 



(frequencies in /iHz). (7) 



Table [T] gives TZ calculated for a set of CoRoT targets with 
solar-like oscillations, with the location f max of the maxi- 
mum PSD and the full-width at half-maximum <5^ onv of the 
pressure mode envelope. The precision of the determination 
of TZ derived from these CoRoT targets is about 15%. 



3.2. EACF as a function of time, filter width and 
signal-to-noise ratio 

The scaling (Eq. [5]) permitted us to perform different treat- 
ments in order to analyze how the EACF varies with the 
observing time t, the filter width Sish, the full- width at half- 
maximum of the mode envelope Su elLV and TZ. With a linear 
dependence of t and the introduction of the reduced width 
X = 5vn/ 5v cm , we found: 



1.5 



t \ 5v em 1 



( * \ 
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0.2 0.4 0.6 0.8 1.0 1.2 1.5 2.0 

Reduced width Su u /6u env 

Fig. 2. Variation of the reduced amplitude 

function of the reduced fil- 
ter width Svu/SveiLv The light grey region indicates the 
validity of the global mean fit given by Eq. ([5]), within 
a ±20% precision. The dark grey region indicates the 
location of the maxima, within a ±15% precision except 
for the double star HD 181906. 



The amplitude Aau increases linearly with time, since a 
is inversely proportional to the observation time t accord- 
ing to Eq. As an important consequence, despite the 
limited lifetime of the modes, the precision of the seismic 
autocorrelation diagnosis increases linearly with t. This in- 
crease corresponds in fact to a decrease of a and there- 
fore cannot saturate. The factor 1.05 in Eq. ([5]) is derived 
from the comparison between a Gaussian envelope and the 
Hanning filter: the best fit with such a filter requires a full 
width at half maximum equal to 1.05 times the one of the 
mode envelope. 

Figure [5] shows the global fit, valid for photometric data of 
solar-like stars obtained with CoRoT or with VIRGO/SPM 
onboard SOHO. All values are fit within ±20%, when Sv-r < 
2Sis env , except the amplitudes for HD175726, which is the 
target with the lowest 1Z\ however, the maximum amplitude 
for this star agrees with the others. 

We have verified that the exponent of the 1Z dependence 
that minimizes the dispersion of the different curves in 
Fig. [2] is 1.5 ± 0.05. A theoretical analysis should be per- 
formed to assess this result. Such a work requires one to 
take into account the link between 1Z and the star inclina- 
tion, the mode lifetime and the stellar noise. 

3.3. Maximum autocorrelation signal 

From Eq. 1)8]) . we can derive the maximum autocorrelation 
signal, obtained for Svn = aSv env . It varies as: 



7.0 a 1Z L 



1 day 



SVr, 



1 mHz 



(9) 



The parameter a is derived from the location of the maxi- 
mum signal (Fig. [2]). If a has not be determined, it should 
be replaced by its typical value a ~ 1.05 as used in Eq. ([8]). 
Figure [5] helps to identify A max . For all solar-like single 
stars but HD181906, the agreement with Eq. © is better 
than ±15%. The fact that HD181906 shows the lowest max- 
imum among solar-like stars is certainly due to its binarity 
(Bruntt 2009): TZ and _4 max are corrupted by the unknown 



for the red giant HD181907, not shown, is 2 times lower 
than expected. This is clearly related to the narrow enve- 
lope of its oscillation spectrum, expressed by 7 = 5 com- 
pared to a mean value of 10 for solar-like stars (Table [1]). 
The number of observed p modes is then twice as small and 
the EACF is reduced. 



4. Performance 

The scaling of *4 max with Eq. © allows us to test the reli- 
ability of the detection of the large separation with the Ho 
test, and then to estimate the scientific output of the EACF. 
The null hypothesis, term fir st coined by th e geneticist and 
statistician Ronald Fisher in Fisher (1935) consists here of 



assuming that the correlation is generated by pure white 
noise. If the EACF is high enough, the Ho hypothesis is 
rejected, implying th at a signal might have been detected 
d Appourchaux 2004[ ). 

4.1. H test 

Assessing the reliability of the measurement of the large 
separation as proposed by RV06 implies applying a sta- 
tistical test as the null hypothesis Ho. A priori informa- 
tion on the large separation may come from scaling laws 
QChristensen-Dalsgaard fc Frandsen 1983[ ) , or may be de- 
rived from the location of the maximum signal, or from 
the initial guess of the stellar fundamental parameters. The 
large separation is then searched for over a range At. The 
number Af of independent bins over the range At depends 
on the width of the cosine filter. It is proportional but not 
equal to the number of points JVh selected by the filter in 
the Fourier spectrum. It can be determined from the full 
width at half-maximum St of the autocorrelation peaks. 
Then, Af is: 



Af 



At 
~5t~' 



(10) 



Therefore the rejection of the Ho hypothesis at probability 
level V implies a threshold value: 



A\\t 



ln(7>) + In 



At 
~St~ 



(11) 



This equation is only valid if V <C 1 . Equation [TT] shows 
that the threshold increases with the searched range, but 
decreases with the resolution St. We verified (sec Eqs. IA.2I 
and IA.3|) that St is related to the width Sv-r of the cosine 
filter. Then, Af is 



Af = - At Svh 



(12) 



with (3 = 0.763. This number Af can be estimated, even 
if nothing is known about the target, since At and 5vn 
are both function of the large separation. The interval At 
where the autocorrelation peak is to be found is measured 
by the time shift t& v . As a conservative value we may con- 
sider At = ta„ = 2/Aza 

As shown above, the width of the best filter giving the max- 
imum autocorrelation signal is proportional to the width of 
the seismic mode envelope, Sua — aSv env , and the mode 
envelope also varies almost linearly with the the large sepa- 
ration, Sv cnv = 7 Av. This gives Svn = «7 Av. As a conse- 
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Table 2. Degree identification 



Fig. 3. Precision of the function Av(v) of HD49933. 1-a 
error bars are given by the dashed lines. The inset shows 
the cosine filter, with a full- width at half-maximum equal 
to 2 times the mean large separation. 




f (rnHz) 

Fig. 4. Same as Fig. El but for HD181420. 

independent bins in the autocorrelation can be estimated 
by: 



T' 



(13) 



Setting the mean optimum value of a to 1.05 and the mean 
value of 7 to 10, we obtain the number of independent bins 
to be considered in the EACF, about 28. The threshold 
values for rejecting the Ho hypothesis at level V — 1% or 
10% are then respectively 8.0 and 5.7. 

4.2. Determination of the mean large separation (Av) 

The determination of the mean value of the large separation 
requires -4 max to be greater than a threshold value of about 
8 for a detection at the 1% rejection level. From Eq. (|A.8|) . 
we then get an estimate of the relative precision of the mean 
large separation (Av), integrated over a large frequency 
centered on the maximum-oscillation frequency: 



8(Av) 



0.057 

^max 



(14) 



At the detection limit -A max = 8, this gives a relative pre- 
cision of about 0.6%. The meaning of a mean value of the 
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T 01 
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HD49385 


7 


2.7 


0.34 


HD49933 IR 


8 


2.0 


2.22 


HD49933 LR 
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6.0 


0.00 


HD181420 


6 


3.1 


0.10 


HD181906 


4 


0.1 


46.6 



IR and LR represent respectiv ely the initial and long runs ob- 
servation on HD 4993 3, from Appourchaux et al. (2008) and 
Benomar et al. (2009)] 



show significantly varying Av(y), These variations, inte- 
grated over the filter width, may limit the precision of (Av) , 
hence of an output value such as the stellar radius. 

4.3. Variation of the large separation with frequency 

With smaller values of 5vn, it is possible to address the 
variation of the large separation with frequency as ex- 
plained in R09. Investigating in detail Av{v) requires a 
filter <5z/h much narrower than <5i^ e nv, so that we can de- 
rive A ~ Amax e/a 5uH./8v em from Eq. ©. In this case, 
the best relative precision at the maximum oscillation fre- 
quency can be derived from Eq. (|A.8|) : 



5Av 
~~AV 



0.6 At/ 
A 5vk 



2.4 /<5i/h 

Anax V Av 



(15) 



The scaling to the large separation insures a uniform pre- 
cision throughout the HR diagram. With <5i/h — 2 Av, 
a 1% precision on the determination of Av(y) requires 
Anax > 60, which is achieved by most of the CoRoT tar- 
gets (Table [J) . Figures [3] and |4] show the precision we can 
obtain on the function Av{y). 



4.4. Disentangling the degree 

Examining the half-separations Aoi = ^,1-^1,0 and A10 = 
^n+1,0 ~ u n,ii a s proposed by R09, requires Svn narrower 
than Av. We have found that 8v-r/ Av = 0.75 provides the 
best compromise: it is narrow enough to select only 1 pair 
of modes with degree and 1, and large enough to give 
an accurate signal-to-ratio. Setting 8vn/ Av = 0.75, a 1% 
precision on the determination of Av requires -4 m ax > 430, 
which is achieved only for HD49933 and HD49385. 
As reported by R09, the different half-large separations are 
clearly distinguished. However, values are correlated within 
the filter, and mixed with other separations including I = 2 
modes. Therefore, we do not consider that the autocorre- 
lation is able to provide a precise measurement of the half- 
separations. For instance, we cannot reproduce the Solar 
values (Fig. However, we clearly show that the local 
minima match the eigenfrcqucncics with high precision. We 
observed, in the unambiguous cases provided by the Sun, 
HD49385 (Fig. [5]) and models, that the local minima asso- 
ciated with 1=1 are lower than the ones with I = 0. This 
is due to the fact that, under the assumption of a Tassoul- 
like spectrum, a narrow filter centered on the 1=1 mode 
tests the separations Aoi and A12 = v n> i — v n ,i, whereas 
a narrow filter centered on the I = mode tests the sepa- 
rations A10 and Aoi- In the asymptotic formalism, A12 is 
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When the filter is not centered on an cigenmode, it mainly 
tests the separation Aoi or Aio- Therefore, we show 
that a narrow frequency windowed autocorrelation allows 
us to distinguish I = from i — 1, which is a cru- 
cial issue since many observations have shown how dif- 
ficult it can be to distinguish them (|Barban c t al. 20091 
IGarcia et al. 2009")) . The test applied to the first initial run 
on HD 49 933 (Fig. [7]) shows that the former mode identi- 
fication of Appourchaux et al. (2008)1 ca nnot be confirmed, 
as also shown by rienomar et al. (2009) who analyze a sec- 
ond longer run. 

A clear identification requires a signal-to-noise ratio high 
enough. Again, Eq. (fT5")) allows us to estimate the auto- 
correlation amplitude required. In order to distinguish the 
small separation D , and considering as a rough estimate 
that in the mean case Dq represents about 2% of the large 
separation, a reliable determination based on a narrow fil- 
ter i5^h = 0.75 A v requires a maximum amplitude greater 
than about 200. 

Table [5] summarizes the mean value of the difference <5oi 
between the local minima compared to the 1-a uncertainty 
SAu/2 of the narrow frequency windowed autocorrelation 
function (Eq. [T5")) : 



An 

E 



°0H 
SAv, 



(16) 



with An the number of pairs of modes above the thresh- 
old level *4 = 8. Table [3] also provides the probability of 
obtaining a result as extreme as the observation assuming 
that the null hypothesis is true. 

This criterion helps to explain why the ridges can be un- 
ambiguously identified in HD 49385 (jDeheuvels et al. 2009|) 
and why scenario 1 for HD 181420 must be preferred 
(jBarban et al. 2009)) : figures [7] and [5] show that the mean 
difference between the local minimal corresponding to 
I = or 1 is greater than the error bar of Av{y). On 
the other hand, no answer can be given for HD 181906 
(|Garcia et al. 2009)) : the low value of A max hampers the 
calculation of Av{v) with a narrow filter (Fig. [9]). The lim- 
ited A max ~ 250 for the initial run on HD49933 and the low 
value 7oi help to explain the difficulties encountered with 
the mode identification given in Appourchaux et al. (2008) 




Fig. 5. Function Av(v) of the Sun, with a filter width equal 
to 0.75 times the mean large separation. The extrema of 
Av(v)/2 do not correspond to the curves Aoi and Aio, 
plotted as dotted curves. The symbols and 1 indicate the 
location of the eigenfrequencies on the A i and Aio curves, 
and indicate also the corresponding local minima of Av(v). 




8v e /<hv> = 0.75 



Fig. 6. Function Av(v) of HD 49385, with a filter width 
equal to 0.75 times the mean large separation. The region in 
grey encompasses the values ± the 1-a error bar. The local 
minima of Av(y), here marked by the degree, correspond to 
the eigenfrequencies mentioned by Deheuvels et al. (2009). 



4.5. Ultimate precision on Av{u) 

Obtaining the best time resolution in the EACF, namely 
the time resolution St of the time series, requires A > 
1.21/&/H St for a filter width Sv^ (cf. Eq. [PO]) . This re- 
lation imposes a strong constraint on the maximum am- 
plitude. Furthermore, in order to investigate the variation 
Av{v), the condition has to be satisfied in a frequency range 
as large as 2Sv cnv around v c . This yields the condition: 



A> 



4.9 
Svn St 



(17) 



According to this, high precision is easier to reach for high 
SiSft values, hence for low mass stars with a higher large sep- 
aration. The decrease of the limit with increasing sampling 



4.6. Small separation 



Roxburgh & Vorontsov (2006) proposed to make use of the 
autocorrelation function to obtain an independent estimate 
of the small separation. This method is based on the com- 
parison of the peak amplitude of even or odd orders in 
the autocorrelation function. With A n the amplitude of the 
peak of order n, it consists of comparing the decreasing A-2n 
values to the increasing Am+\ (see Fig. 4 of RV06). Note 
that Aav corresponds to A2 ■ Equality of the interpolated 
curves A2n and A2 n +i occurs for n of 3 or 4. Tests made on 
the available C0R0T data show that values of Ain larger 
than 15 are necessary to apply the method, therefore re- 
quiring very large values of Aau, larger than about 300. 
In the solar case, with the simulation including additional 
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HD -19933 IR 




Su B = 63 /zHz 
Sv n /<Av> = 0.76 




Fig. 7. Function Au(u) of HD49933, with a filter width 
equal to 0.75 times the mean large separation. Radial and 

1 = 1 modes are identified. The two plots correspond to the 

2 runs: the initial run IR lasted 60.7 days and the long run 
LR 136.9 days. 
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6l> H = 55 /^Hz 
Sv K /<At/> = 0.74 
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Fig. 8. Same as Fig.0 but for HD181420. 



4. 7. Threshold levels 

Table [3] summarizes the threshold levels for the determi- 
nation of the seismic parameters with the EACF. At low 
signal-to-noise ratios, namely an A max value lower than 8, 
the method cannot operate, according to the Ho test. Then, 
the domain where the autocorrelation is highly perform- 
ing is for A max ranging from 8 (detection limit) to ~ 50, 
when precise mode fitting becomes possible (HD181906, 



3. 

« 42 



HD181906 



if 




5v H = 83 /iHz 
6v B /<&!/> = 1.00 



1.6 1.7 l.B 1.9 2.0 

f (mHz) 

Fig. 9. Same as Fig. [3 but for HD181906 and with a 
broader filter. The large uncertainty indicated by the broad 
grey region shows that the identification for that star is not 
reliable. 

Table 3. Threshold levels 



An 



detection 



< 5.7 no detection (10% rejection level) 

< 8.0 no detection (1 % rejection level) 

10 measurement of (Ai/) with 0.5 % precision 
> 50 fitting of the modes has proven to be possible 

> 200 identifying the mode degree has proven to be pos- 

sible 

> 300 possible estimate of the smal l separation as in 
Roxburgh fc Vorontsov (2006) | 



The threshold values are given for a typical F dwarf. They can be 
made precise for a given star, according to its specific parameters 
a and 7. Levels are lower for red giants (to be estimated with 
7 ~ 4 instead of 10). 



HD49933 
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3 6 12 24 48 96 192 

2/T 4li ( M Hz) 

Fig. 10. Automatic search for the signature of a large sepa- 
ration for HD49933. The grey line indicates the location of 
the maximum peak, corresponding to the mean large sep- 
aration of the star. The horizontal segments indicate the 
ranges corresponding to the 13 initial guess values Au c . 



may be useful for identifying the degree of the modes, un- 
der the condition that the oscillation spectrum is close to a 
Tassoul-likc pattern. Larger A max values allow a more de- 
tailed analysis with classical methods such as mode fitting 
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Table 4. Kepler performance 
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for Kepler performance on CoRoT-like targets with varying magnitudes. 
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Fig. 11. Same as Fig. [10l but for HD175726. The dashed 
line indicates the 10% rejection limit. 
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Fig. 12. Same as Fig. El but for the red giant HD181907. 
The horizontal dashed lines indicate the 10% and 1% re- 
jection limits. The vertical grey lines indicate the signature 
at Av and 2 Av, the dotted line the spurious signatures of 
the day aliases, and the dot-dashed lines the CoRoT orbital 
frequency and first subharmonic. 



5. Discussion 

5.1. Automated determination of the large separation 

Autocorrelation may provide an effective automatic deter- 
mination of the large separation when nothing is known 
about the star, as can be the case for a Kepler target. As 
shown previously, testing the autocorrelation around ta^ 
requires a cosine filter Svn = o-l Av ~ 10.5 Av, near the 



conditions, in a frequency range simultaneously encompass- 
ing giant and dwarf stars, i/ ma j is derived from the scaling 
law given by Eq. ([7]). 

The automatic test consists of analyzing the autocorre- 
lation of the time series for a set of time shifts t& vc in 
geometrical progression. We performed the automatic au- 
tocorrelation test with 13 values of Av c = 2/tau c , vary- 
ing from 3 to 192 /iHz with a geometric ratio G equal to 
\/2; Av c can be considered as an initial guess of the large 
separation. For each initial value Av c , we explored the 
range [Av c /G, Av c G] of the autocorrelation for 3 frequency 
ranges of the Fourier spectrum centered respectively around 
i/ mal and iAnax i 5isn/2. We finally derived the large sep- 
aration from the maximum amplitude .4 au to calculated for 
each Av c initial guess. Comparison of the different -4auto is 
made possible by the scaling provided by Eq. ©. Figure ITTJ1 
shows the result for HD49933. 

We also tested the automatic test with the stars with 
the lowest A max , namely HD175726 and HD181907. For 
HD175726, the single value exceeding the 10% rejection 
level occurs at 97 /iHz (Fig. QT]). This value of Av agrees 
with the solution proposed by Mosser et al. (2009). This de- 
tection is poor since a significance level of 10 % means that 
the posterior proba bility of the null hypothesis is at least 
38 % according to Appourchaux et al. (2009) However, 
the automatic detection can be refined with a dedicated 
search with a more precise grid of analysis. In the case of 
HD175726, the clear identification of an excess power cen- 
tered at 2 mHz first allows us to better estimate the parame- 
ters for searching Av and, second, gives a further indication 
that the measurement is reliable thanks to Eq. (JT)). 
The amplitude -4 au to of the automatic test is found to be 
close to the maximum amplitude *4 max . Only limited fine 
tuning around the automatically fixed parameters is needed 
to optimize the result. Mosser et al. (2009) have mentioned 
the difficulty of determining the large separation with other 
methods. The autocorrelation method proves to be pow- 
erful for a rapid estimate of the large separation; rapidity 
means here a few seconds of CPU time for the Fourier spec- 
trum of the CoRoT time series, followed by a few seconds 
of CPU time for the automatic search with the autocorre- 
lation, with a common laptop. 

We verified that the method is effective for all other solar- 
like targets: it gives one single answer for Av, and does not 
deliver any false positives. The case of red giants requires 
a dedicated analysis. 

5.2. Red giants 

We tested the method on the CoRoT red giant target HD 
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Fig. 13. Large separation, auto matically measured for a 
set of 392 red giants analyzed in Hekker et al. (2009)[ as a 
function of the maximum oscillation frequency i/ max . Large 
black squares indicate positive detection; small grey squares 
correspond to unreliable cases. The dashed line indicate the 
global fit described by Eq. [TH a few cases correspond to the 
automatic detection of twice the large separation. 



get (about 2) and a large A max autocorrelation signal, the 
large separation is easily found around 3.5 //Hz, but the 
detection is polluted by many values clearly above the 1 % 
rejection limit (Fig. [T2"|) . All these spurious detections are 
caused by artefacts: detection of the double of the large 
separation; detection of the diurnal frequency and its har- 
monics; detection of the CoRoT orbital frequency and half 
its value. We checked that the detections at high harmonics 
of the diu rnal frequency are d ue to residuals of the window 
function (M osser et al. 2009)) ; they are introduced by the 
link between f max and Av indicated by Eq. [7] 
For HD 181907, Eq. © is valid within 30%. As discussed 
in Section [31 the discrepancy compared to solar-like stars 
is due to the fact that the mode envelope of red giants is 
narrower than in solar- like stars. 

The automated determination of the large frequency has 
been also tested on a set of 392 giants observed in the 
CoRoT fi eld dedicated to exo planetary science and ana- 
lyzed by Hekker et al. (2009)| The method proves to be 
efficient and rapid. It provides a clear advantage since it 
gives a quantified reliability thanks to the use of the Ho 
test. We present in Figure [TBI the results obtained for these 
gi ants. The maximum o scillation frequency was calculated 
by Hekker et al. (2009) The amplitude of the autocorrela- 
tion signal allows us to clearly discriminate artefacts from 
reliable detection (60% of the targets). The relative preci- 
sion in the mean large separation is much better than 1 %, 
according to Eq. [I4j After correction of the stars for which 
the double of the large separation is preferably automati- 
cally detected, we define from the fit of the relation between 
the large separation and the location of the maximum sig- 
nal a power law varying as: 



5.3. Kepler data 

The Kepler mission compared to CoRoT will provide 
different photometric performance, on dimmer targets 
but in some cases with longer observation duration 
( jChristensen-D alsgaar cTet al. 2007| . According to Kepler 
performance (Kjeldsen et al. 2008), the noise level is about 
0.92, 10.2 and 144ppm z /iHz _i for targets of V magnitude 
respectively equal to 9, 11.5 and 14. 

We can extrapolate the performance obtained with Kepler 
on targets similar to the ones observed by CoRoT, but of 
magnitude 9 to 14, after a 4- year long observation. Tabled 
gives the amplitude -4 max for targets observed during typ- 
ical 90-day or 4-year long runs. According to the expected 
performance in 90-day runs, the brightest F-type or the 
class IV targets will have a signal-to-noisc ratio high enough 
to derive information on the large separation. In a 4-year 
run, the brightest G dwarfs will deliver a clean seismic sig- 
nature. On the other hand, faint F targets will have fully 
exploitable Fourier spectra that will require a precise mode 
fitting for the most complete seismic analysis. The perfor- 
mance for giants appears to be almost independent of the 
magnitude, since the contribution of photon noise is negli- 
gible at low frequency. 

We can compare this approach to the hare -and-hounds ex- 
ercises performed by Chaplin et al. (2008) The asteroseis- 
mic goal of Kepler is principally to derive information on 
stars hosting a planet, by the determination of the large 
separation. Compared to global fitting, the autocorrelation 
function gives a more rapid and direct answer. 
The autocorrelation benefits from the rapid cadence (32 
s) provided by CoRoT in the seismology field. Kepler will 
provide 2 cadences, at 1 or 30 min. This yields a lower res- 
olution in time, hence a lower precision on the expected 
results. 



6. Conclusion 



Roxburgh & Vorontsov (2006) have proposed a method for 
estimating large and small sepa rations from the analysis 
of the autocorrelation function. Roxburgh (2009) has ex- 
tended the method to determine the variation of the large 
separation. In this paper, we have developed and quantified 
the method, relating the amplitude of the correlation peak 
at time shift t^ v = 2/ Av to various parameters. 
We have scaled the autocorrelation to the white noise con- 
tribution, so that we were able to relate the autocorrelation 
signal to the mean seismic height-to-background ratio 1Z 
that measures the relative power density of the signal com- 
pared to noise and to background signals. This empirical 
relation is precise to about 15% for solar-like stars. 1Z ag- 
gregates the influence of unknown parameters such as the 
mode lifetimes, the star inclination (that governs the modes 
visibility) or the rotational splitting. On the other hand, all 
these unknown parameters complicate and slow down the 
fitting of individual eigenfrequencies. Therefore, the EACF 
shows here a possible advantage in terms of speed. 
The EACF gives a direct measurement of the mean large 
separation. Compared to other methods, the estimate is ac- 
curate and simple, with an intrinsic threshold value, with 
error bars, and without any modeling of the other com- 
ponents of the Fourier spectrum (granulation or activity). 
This law for giants is in agreement with Eq. (IT]) based on Furthermore, when the signal-to-noise ratio is high enough, 



Av~ (0.26 ±0.015) v, 



0.78±0.03 



(frequencies in /iHz). (18) 
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large separation with frequency, without any mode fitting. 
This is a key point for stellar radius measurement. Previous 
works have shown the difficulty to disentangle I = from 
I = 1 modes in oscillation spectra of F stars observed in 
photometry. We have verified that, for high signal-to- noise 
ratio Fourier spectra, the autocorrelation analysis can pro- 
vide an unambiguous identification of the mode degree for 
a solar-like oscillation spectrum. 

We have defined a method for the automatic determination 
of the large separation, which is efficient at low signal-to- 
noise ratio, even if no information is known for the star. 
We have determined that the width of the cosine filter used 
in the method that optimizes the EACF is very close to 
full-width at half-maximum of the mode envelope (ratio 
about 1.05). We have also checked that the performance 
of the method increases linearly with the duration of the 
time series. With very limited CPU time (a few seconds), 
this method delivers the mean large separation of a target. 
It requires no information on the star; it just relies on the 
assumption that the location of the excess power and its 
width are related to the large separation by a scaling law, 
what is verified for red giants and solar-like stars. Finally, 
we were able to investigate in a simple manner the capabil- 
ity of Kepler. 

We are confident that the autocorrelation method will 
be of great help in analyzing high duty cycle time series 
as a complement to the Fourier analysis. As noticed by 
Fossat et al. (1999)] the autocorrelation signal gives a clear 
signature since the autocorrelation delay, namely four times 
the stellar acoustic radius (about 4 to 8 hours for an F 
dwarf), is much shorter than the mode lifetime (a few days). 
This allows each wavepacket to properly correlate with it- 
self after a double travel along the stellar diameter, so that 
the autocorrelation integrates phased responses over the 
total duration time. In the Fourier spectrum, on the con- 
trary, interference between the short-lived wavepackets ob- 
served in the time series produce a complicated pattern. 
But Fourier analysis still remains required for the precise 
determination of the eigenfrequencies derived from an ac- 
curate mode fitting. 

Appendix A: Performance of the autocorrelation 

A.l. Square module of the autocorrelation 

The EACF presented in Section [2] is defined to directly 
give the envelope of the autocorrelation. Since negative fre- 
quencies are omitted, the EACF is related to the canonical 
autocorrelation C± , that includes positive and negative fre- 
quencies of the Fourier spectrum, by: 



Ak\C\* and \C\ = JC± Z +H(C±) 



S 0.5 



o.o - 



(A.l) 



where ft is the Hilbert transform. Fig. IA.1I shows the dif- 
ference between C and C± , both functions include the con- 
tribution of a narrow Hanning filter. 

A. 2. Autocorrelation peak 

The shape of the autocorrelation peaks is given by 
the Fourier transform of the Hanning filter, which 
can be expressed a s the sum of 3 components 
(|Max fc Lacoume 1996)) : 




-0.5 



Fig. A.l. Comparison of ± |C| (black curves) and C± (grey 
line), both normalized to 1 at time shift 0, calculated for 
the narrow frequency windowed solar Virgo spectrum. 



with sincA = sin(7rA)/7rA. This global shape may be more 
simply modeled with a single cosine shape (Fig. IA.2[) . In 
order to enhance the precision of the fit in the upper part 
of the autocorrelation peak, the full width at half maximum 
St of this cosine fit is: 



5r = -f- with 13 ~ 0.763. 



(A.3) 



The flanks of the peak are not well fitted, which is unimpor- 
tant compared to the fact that the fit above half-maximum 
performs well. St is much greater than the resolution time 
St. 

Precise determination of the large separation requires pre- 
cise location of the peak maximum. In order to estimate 
the performance, we describe a peak as: 




1 



COS 7T- 



_* — TAu 

St 



for \t\ < St, 
for \t\ > St. 



(A.4) 



This fit shows variation 
dS A 



TtA t- T Au 

—— sm7r dr. 

2St St 



(A.5) 



We can compare the variation of the signal peaking at am- 
plitude A to the maximum variation of a noise contribution 
of amplitude b. At a time shift At from the maximum, the 
signal variation and the maximum noise contribution are: 



AS* 



AS h 



A 
2 
b 

2 



At 

~~Sr~ 
At 

J?' 



(A.6) 



The precise identification of the signal maximum requires 
A«Sa > ASb, which translates into the condition: 



7T A At > b St. 



(A.l) 
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Fig. A. 2. Autocorrelation peak A/\ v * (solid line) normal- 
ized to 1 at time shift 0, theoretical profile (dashed line) and 
cosine model (dotted line) calculated from the HD 49933 
power spectrum and a filter width Sv-r = 0.2 mHz. 

A. 3. Precision on the mean large separation 



If the amplitude is not large enough, then Eq. IA.7I defines 
a resolution At, hence a limited precision on the large sep- 
aration: 



SAv 
~Ai7 



At 

TAv 



(3 b Av 
2?r A Svh 



(A, 



We can set, as a limit to detection, 6=5 and A = 
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A. 4. Full resolution for the measurement of Av(v) 

In order to recover the full time resolution (At = St), the 
amplitude must satisfy A > Ast, with the definition: 

As - -— - hf3 1 (A 9) 

TT St TT Svft St 

With 6 = 5: 

Ast * T~~T7 • (A.10) 
Sun St 
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